21 research outputs found
Sapporo2: A versatile direct -body library
Astrophysical direct -body methods have been one of the first production
algorithms to be implemented using NVIDIA's CUDA architecture. Now, almost
seven years later, the GPU is the most used accelerator device in astronomy for
simulating stellar systems. In this paper we present the implementation of the
Sapporo2 -body library, which allows researchers to use the GPU for -body
simulations with little to no effort. The first version, released five years
ago, is actively used, but lacks advanced features and versatility in numerical
precision and support for higher order integrators. In this updated version we
have rebuilt the code from scratch and added support for OpenCL,
multi-precision and higher order integrators. We show how to tune these codes
for different GPU architectures and present how to continue utilizing the GPU
optimal even when only a small number of particles () is integrated.
This careful tuning allows Sapporo2 to be faster than Sapporo1 even with the
added options and double precision data loads. The code runs on a range of
NVIDIA and AMD GPUs in single and double precision accuracy. With the addition
of OpenCL support the library is also able to run on CPUs and other
accelerators that support OpenCL.Comment: 15 pages, 7 figures. Accepted for publication in Computational
Astrophysics and Cosmolog
The dynamics of stellar disks in live dark-matter halos
Recent developments in computer hardware and software enable researchers to
simulate the self-gravitating evolution of galaxies at a resolution comparable
to the actual number of stars. Here we present the results of a series of such
simulations. We performed -body simulations of disk galaxies with between
100 and 500 million particles over a wide range of initial conditions. Our
calculations include a live bulge, disk, and dark matter halo, each of which is
represented by self-gravitating particles in the -body code. The simulations
are performed using the gravitational -body tree-code Bonsai running on the
Piz Daint supercomputer. We find that the time scale over which the bar forms
increases exponentially with decreasing disk-mass fraction and that the bar
formation epoch exceeds a Hubble time when the disk-mass fraction is
. These results can be explained with the swing-amplification theory.
The condition for the formation of spirals is consistent with that for
the formation of the bar, which is also an phenomenon. We further argue
that the non-barred grand-design spiral galaxies are transitional, and that
they evolve to barred galaxies on a dynamical timescale. We also confirm that
the disk-mass fraction and shear rate are important parameters for the
morphology of disk galaxies. The former affects the number of spiral arms and
the bar formation epoch, and the latter determines the pitch angle of the
spiral arms.Comment: 23 pages; 29 figures. Accepted by MNRA
A distributed SIRT implementation for the ASTRA Toolbox
The ASTRA Toolbox is a software toolbox that enables rapid development of GPU accelerated
tomography algorithms. It contains GPU implementations of forward and backprojection operations
for common scanning geometries, as well as a set of algorithms for iterative reconstruction. These
algorithms are currently limited to using a single GPU.
A drawback of iterative
reconstruction algorithms is that they are slow compared to classical
backprojection algorithms. As a result, using only a single GPU can result
in prohibitively long reconstruction times when working with large data volumes.
In this paper, we present an extension of the ASTRA Toolbox with implementations of forward
projection, backprojection and the SIRT algorithm that can be distributed over
multiple GPUs and multiple workstations to make processing larger data sets
with ASTRA feasible
A distributed ASTRA toolbox
While iterative reconstruction algorithms for tomography have several advantages compared to standard backprojection methods, the adoption of such algorithms in large-scale imaging facilities is still limited,
Trimodal structure of Hercules stream explained by originating from bar resonances
Gaia Data Release 2 revealed detailed structures of nearby stars in phase
space. These include the Hercules stream, whose origin is still debated. Most
of the previous numerical studies conjectured that the observed structures
originate from orbits in resonance with the bar, based on static potential
models for the Milky Way. We, in contrast, approach the problem via a
self-consistent, dynamic, and morphologically well-resolved model, namely a
full -body simulation of the Milky Way. Our simulation comprises about 5.1
billion particles in the galactic stellar bulge, bar, disk, and dark-matter
halo and is evolved to 10 Gyr. Our model's disk component is composed of 200
million particles, and its simulation snapshots are stored every 10 Myr,
enabling us to resolve and classify resonant orbits of representative samples
of stars. After choosing the Sun's position in the simulation, we compare the
distribution of stars in its neighborhood with Gaia's astrometric data, thereby
establishing the role of identified resonantly trapped stars in the formation
of Hercules-like structures. From our orbital spectral-analysis we identify
multiple, especially higher order resonances. Our results suggest that the
Hercules stream is dominated by the 4:1 and 5:1 outer Lindblad and corotation
resonances. In total, this yields a trimodal structure of the Hercules stream.
From the relation between resonances and ridges in phase space, our model
favored a slow pattern speed of the Milky-Way bar (40--45 ).Comment: 11 pages, 9 figures, MNRAS accepte
Impact of bar resonances in the velocity-space distribution of the solar neighbourhood stars in a self-consistent -body Galactic disc simulation
The velocity-space distribution of the solar neighbourhood stars shows
complex substructures. Most of the previous studies use static potentials to
investigate their origins. Instead we use a self-consistent -body model of
the Milky Way, whose potential is asymmetric and evolves with time. In this
paper, we quantitatively evaluate the similarities of the velocity-space
distributions in the -body model and that of the solar neighbourhood, using
Kullback-Leibler divergence (KLD). The KLD analysis shows the time evolution
and spatial variation of the velocity-space distribution. The KLD fluctuates
with time, which indicates the velocity-space distribution at a fixed position
is not always similar to that of the solar neighbourhood. Some positions show
velocity-space distributions with small KLDs (high similarities) more
frequently than others. One of them locates at , where and are the distance from the galactic centre
and the angle with respect to the bar's major axis, respectively. The detection
frequency is higher in the inter-arm regions than in the arm regions. In the
velocity maps with small KLDs, we identify the velocity-space substructures,
which consist of particles trapped in bar resonances. The bar resonances have
significant impact on the stellar velocity-space distribution even though the
galactic potential is not static.Comment: 9 pages, 11 figures. Accepted by MNRA